Mon. Not. R. Asti'on. Soc. 000,[JJ-?? (201 1) Printed 18 January 2012 (MN KTgX style file v2.2) 



<n : 



o 

u 

Oh 
6 

o3 



> 

m 
cn 

o 

(N 



On the dust abundance gradients in late-type galaxies: 

II. Analytical models as evidence for massive interstellar dust 

growth in SINGS galaxies 

Lars Mattsson* & Anja C. Andersen 

Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100, Copenhagen 0, Denmark 



ABSTRACT 

We use simple analytical models of the build up of the dust component and compare these with 
radial dust distributions derived from observations of SINGS galaxies. The observations show 
that dust gradients are indeed typically steeper than the corresponding metallicity gradients 
and our models indicate very little dust destruction, but significant dust growth in the ISM for 
most of these galaxies. Hence, we conclude that there is evidence for significant non-stellar 
dust production, and little evidence for dust destruction due to SNe shock waves. We find 
that dust is reprocessed rather than destroyed by shocks from SNe. Finally, we argue that dust 
abundances derived using standard methods may be overestimated, since even very 'generous' 
estimates of the metallicity results in dust-to-metals ratios above unity in several cases, if the 
dust abundances given in the literature are taken at face value. 
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1 INTRODUCTION 

Dust grains in dense molecular can theoretically grow to microm- 
eter sizes by accre t ion and coagulation (see , e.g.. IOssenkopl][r993l : 
lOrmel et all 120991 : Iffirashita & Kuol l201lh . There is also obser- 
vational evidence suggesting such l arge grains are abu ndant in 
many Galactic molecular clouds (see lPagani et alJ[2010L and ref- 
erences therein), which is easily explained by efficient dust growth 
in these cores. Micrometer-sized dust grains can theoretically be 
forme d in carbon rich AGB stars as well (Mattsson & Hofnei 
1201 ll). and ar e likely to form also in oxygen rich AGB stars 
faofned 12008). but the probability of having any larger quanti- 
ties of these grains gathering in dense molecular cores is very low. 
Furthermore, grains are expected to survive in the Galactic ISM 
for a few hundred Mvr {jones. Tielens & Hollen bach 1996; Jones 
|2004 ISerra Dfaz-Cano & Jonesll2008l ; Ijones & NutbH201 ll) . while 



the time scale for any significant stellar dust enrichment in the local 
ISM is about Gyr, which implies the presence of dust growth pro- 
cesses in the ISM. Hence, interstellar dust growth can be regarded 
an established phenomenon. 

Shock-waves from supernovae (SNe) is thought to be able 
to destroy dust grains as these waves propagate through the in- 
terstellar medium (ISM). The time scale for this dust destruction 
is largely dependent on the supe rnova rate (SNR) and how much 
dust each SN-shock can destroy dMcKedl 1 9891: iDrainell 1 990h . This 
shock destruction of dust grains has been conside red quite efficient 
jjones. Tiele ns & Hollenbach 1 19961 : Ijonej [2004), although more 
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recent work show that while this i s likely the case for carbon dust 
it may not be the cas e for silicates JSerra Dfaz-Cano & Jonesl 2008: 
|jones&Nuthll201lh which seems consistent with the Milky Way 
(solar neighbourhood). Very efficient dust destruction seems on 
the other hand i nconsistent with the high dust masses detected in 
high-z objects ( lDweketal.ll2007l : lGall, Andersen & Hiorthll201ll: 
Mattsson 2011). It may of course be that dust destruction by SNe is 
less efficient in high-z galaxies, but also modelling of nearby late- 
type galaxie s seems to work nicely witho ut significant net destruc- 
tion of dust (Inoue 2003; H irashi tall 19991) . and it makes it easier to 
explain the dust-to-gas ratios, since stellar dust production is insuf- 
ficient fairashitdll999l : IZhukovska. Gail & Trieloffll2008h . 

Obser vations imply that dust p roduction in SNe is rather 
inefficient dKotak et al.l l2006t | 2009), which suggest the large 
masses of cold dust detected in some SN remnants (see, e.g. , 
iMorgan & Edmunddl200l : iMorgan et al]|2003l ; iDunne et alj|2009h 
may be the result of subsequent dust growth rather than an 
effect of heating or actual dust production in the actual SN. 
Theoretical results suggest 90% of the dust produced in SNe 
is destroyed by the rever se shock before it reaches the ISM 
jBianchi & S chneider 2007), which leads to a relatively consistent 
pictur e where AGB stars is the most important sour ce of stellar 
dust jEdmunds & Ealesll 19981 : iFerrarotti & Gailll2006f). Mode ls of 
the dust evolution in the solar neig hbourhood bv lDwekl l l 19981) and 
IZhukovska. Gail & Trielofij d2008l) suggest dust production in stars 
and possible dust destruction by SN-shock waves need to be com- 
pensated by an efficient dust growth in the ISM. This is also sup- 
ported by the fa ct that dust grains of considerable sizes seem to 
exist in the ISM (Pagani et al. 2010). However, recent results from 



©2011 RAS 



2 Mattsson & Andersen 



Herchel observations of SN 1987A show a huge reservoir of cold 
dust which sugg est that at least some S Ne may produce significant 
amounts of dust (Matsuura et alj|20llh . However, it is still unclear 
how much of this dust will survive once the reverse shock sets in. 

It is obviously not clear whether dust destruction from SN- 
shock waves in the ISM is an efficient mechanism or not, or per- 
haps more accurately put: whether the net dust destruction is im- 
portant or not. As mentioned above, existing models of the evolu- 
tion of the dust phase in the solar neighbourhood (solar circle) rely 
on the assumption that both growth and destruc tion are important 
dDweklll998l ; IZhukovska, Gail & Trielofll koOS). But with only a 
single data point, e.g., the present-day dust abundance in the solar 
neighbourhood, it is hard to tell whether a model with, or with- 
out, dust destruction by SNe and/or dust growth in the ISM makes 
a better fit to data. Recently , however, results from the SINGS 
project dKennicutt et alj|2003l) have provided dust abundance pro- 
files base d on physical dust models for a good number of nearby 
galaxies dMunoz-Mateos et ai1l2009l) and the associated THINGS 
proj ect has provided H i profiles for many of the SINGS galax- 
ies dWalter et alj|2008h . Hence, the dust-to-gas ratio along the disc 
can be derived for many of these galaxies, which opens for more 
precise comparisons with observations since different galactocen- 
tric distances in the discs represents different evolutionary states of 
the interstellar gas. Many of the late-type galaxies in this sample 
have also well-constrained metallicity gradients as there are suf- 
ficient numbers of detectable Hn-regio ns for which spectroscopy 
have been done (Moustak as et al]20ld) . 

Previous work on dust evolution modelling has shown that the 
dust-to-metals ratio may change as the evolut i on pro ceed s. The an- 
alytica l relations derived by, e.g., Edmunds (2001) and iMattssonl 
J20I 1), show that the dust abundance may not necessarily fol- 
low the metal abundances in the ISM. More detaile d results, ob- 
tained through numerical modelling, suggest th e same dDwekl 19981 ; 
Ilnoudl2003l:lzhukovska. Gail & Trielofl12008l) . The reasons for the 
changing dust-to-metals ratio are mainly the destruction (including 
astration) and growth of dust grains that may occur in the ISM, 
which suggest that dust and metallicity gradients may not have the 
same slope, i.e., there is a dust-to-metals gra dient. As suggested in 
iMattsson. Andersen & Mun khammar l i201ll . henceforth 'Paper F) 
the dust-to-metals gradient along a galactic disc can essentially be 
regarded as a diagnostic for net dust growth or net destruction de- 
pending on whether the slope is negative or positive, respectively. 

We demonstrate here a way to estimate the efficiency of the 
net dust destruction by SNe in nearby late-type galaxies using dust- 
and metallicity gradients. Comparing the dust abundance gradients 
along the disc with the metallicity gradients may give an indication 
of whether net destruction or net growth of dust is important or not, 
as we will show in this paper. 



2 MODEL 

2.1 Basic equations 

Following Paper I we use the instantaneous recycling approxima- 
tion (IRA, which essentially means all stars are assumed to have 
negligible lifetimes compared to the time scale for the build-up of 
metals and dust, see IPagelll 19971) throughout this paper. Then, as- 
suming a 'closed box', gas and dust destruction in the ISM from 
SN-shocks, the equations for the metallicity Z and the dust-to-gas 
ratio Zd becomes 
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where Td is the dust destruction time scale, r gr is the dust-growth 
time scale and y$, y z denote the yields (dust and metals, respec- 
tively) as defined in Paper I, and S s , S s denotes surface density by 
mass of stars and gas, respectively. Combining equations l[T) and 
l(2j, we have 



dZ d _ Jd +Z d (r gr 1 -T d ') 
dZ y z 



(3) 



which thus have no explicit dependence on the gas mass density E e 
or the stellar mass density S s , although Z, Zj as well as the time 
scales T 5r , r d may implicitly depend on X g and S s .. 

2.2 Destruction and growth of dust in the ISM 

Following lMcKed dl989l) the dust destruction time-scale is 



(4) 
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where (otism) is the effective gas mass cleared of dust by each SN 
event (which is not necessarily the same as the gas mass affected 
by each SN), and R$ N is the SN rate. As shown in Paper I, the time 
scale Tj may be expressed as 

-1 „ S d ^ 



(5) 



where S will be referred to as the dust destruction parame- 
ter, which is a kind of measur e of the efficiency of dust de- 
struction. For a lLarsor] d 19981) I MF and to B m « 1000M o 
dJones. Tielens & Hollenbachll 19961: |jonesll2004l) . then d ~ 10 (see 
iMattssorfeOllh . 

We define the rate per unit volume at which the number of 
atoms in dust grains m ay grow by a ccretion of metals onto these 
dust grains as (see, e.g. lDweklll998l) at a rate which naturally de- 
pends on metallicity Z, the typical grain radius a and the stick- 
ing coefficient / s (i.e., the probability that an atom will stick to the 
grain). Additionally, the mean thermal speed of the gas particles 
(including metals) (v g ) will affect this rate, which suggest a tem- 
perature dependence. The timescale of grain growth can thus be 
expressed as (again, see Paper I for brief derivation) 
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where (m gr ) is the average mass of individual dust grains in the 
gas, d c is the characteristic size of the molecular clouds where 
dust grows and S mol is the surface density of molecular gas. As 
in Paper I, we will assume S mo i « Eh 2 , since most of the gas in 
molecular gas clouds is in the form of molecul ar hydrogen. We 



also as sume £ H o trace s the s ta r-formation r ate (Rownd & Young 
1 19991 : IWong&Britzl 120021: IBigiel et alj l2008l: iLerov eTaT 



20081: IBigiel et al.1 bOl il TfFeldrnann. Gnedin & Kravtsovl I20T1F 
ISchruba et al.ll201 11) . This is als o supported by theory and recent 
nume r ical experiments (see, e.g. | Krumholz. McKee & Tumlinsonl 



2009; iKrumholz, Lerov & McKed l201ll) . Moreover, the mean 
thermal speed (v g ) is roughly constant in the considered ISM 
environment and the typical grain radius does not vary much. 
Hence, as shown in Paper I, To is essentially just a simple function 
of the metallicity and the growth rate of the stellar component, 
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S„ dt 



(7) 



where e will, in the following, be treated as a free (but not uncon- 
strained - see Paper I) parameter of the model. The rate of change of 
the dust-to-gas ratio Z d due to accretion of metals onto pre-existing 
dust grains in the ISM is then 



\dt L \ Z)t p 



(8) 



The e used above should also not be confused wit h the dust- 
desruction efficiency parameter used by |McKed ( fl989T) and it also 
differs slightly from the dust-growth parameter used in iMattssonl 
J20l 1) as we in this work consider the free atomic metals rather 
than the total metallicity. 



2.3 Models of radial dust-to-gas distributions 

For simplicity we stick to a closed-box model as in Paper I, i.e., 
no in- or outflows to/from the disc will be considered. This is not 
in agreement with the widely accepted ideas about galaxy-disc for- 
mation, where the baryons (in the form of essentially pristine gas) 
are assumed to be accreted o ver an extended period of time. But 
as shown bv lEdmundsl dl99Cj|) . the most prominent effect of unen- 
riched infall is to make the effective yield smaller, i.e., to dilute the 
gas so that the metallicity builds up more slowly. As we use the 
present-day metallicity as input (see section [3"72l . the overal l effects 
of ass uming a closed box are rather small. Note also, that lGarnettl 
d 19971) have showed that the simple closed-box model of chemi- 
cal evolution seem to work quite well for modelling oxygen (O/H) 
gradients. In particular, adopting the Clavton ( 1987) infall model 
neither means a significant improvement of the fit, nor a very large 
change of the required yield. 

Adopting a closed-box scenario, the dust destruction and dust 
growth models as described above, we arrive at the equation 



§^(-4(-f)M)- 



(9) 



where y z is the metal yield. With < y<j < y z as a basic require- 
ment, the general closed-box solution (of equation [9]l for the dust- 
to-gas ratio Zd in terms of the metallicity Z can then be expressed in 
terms of the confluent hypergeometric Kummer-Tricomi functions 
of the first and second kind (usu ally denoted U and M), respectively 
jKummeJl 1 8371 : ITricomi 1947). We refer to Paper I for further de- 
tails. In case there is no dust destruction by SNe (S = 0) the solution 
reduces to 
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where M(a, b;z) = {F\ (a, b; z) is the Kummer-Tricomi function of 
the first kind, which is identical to the confluent hypergeometric 
function iFi (a, b; z) that is often implemented in computer algebra 
software. If there is neither growth, nor destruction of dust in the 
ISM (e = S = 0), we have the trivial case 



y -±z, 

yz 



(ii) 



corresponding to pure stellar dust production. These are the two 
cases we consider in detail in this paper. 



3 RESULTS AND DISCUSSION 

We have compared the theory of interstellar dust evolution pre- 
sented in Paper I with data on dust gradients in late-type galaxies 
taken from the analysis of SINGS sample bv lMunoz-Mateos et al.l 
(2009). They present radial distributions of dust (surface densi- 
ties) in galactic discs derived from UV-, IR- and gas mass profiles, 
which have relatively small errors. These data are likely the best 
constraints on radial dust gradients in galactic discs available at 
present. Here we summarise and discuss the results. 



3.1 General results and limitations 

In a previous paper (Paper I) we have shown that dust destruction 
by shocks from exploding SNe will flatten a negative (declining 
with galactocentric distance) dust-to-gas gradient over a galactic 
disc, while dust growth in general acts as to make it steeper (see 
the Theorem proved in Paper I). From this result, we expect dust- 
to-metals gradients to have positive gradients (rising with galac- 
tocentric distance) if dust destruction is more important than dust 
growth, and if dust growth is the more important process we expect 
them to be negative. 

In our limited samp le of late-type galaxies taken from 
iMunoz-Mateos et al. we have found that the dust-to-metals 

gradients are typically negative or in some cases flat. There is only 
one counter-example (NGC 5194) where the metallicity gradient 
may be somewhat steeper than the dust-to-gas gradient (see fig- 
ure Q3, resulting in a positive dust-to-metals gradient (see figure 
[2). However, it should be stressed that the metallicity gradients (as 
traced by the O/H gradients in the discs) we adopt in this study 
may not be correct in all cases. We have taken the metallicity gra- 
dients from Moustak as et al.l d201dh. which are derive d using the 
strong-line calibration bv lKobulnickv & Kew lev (2004), except for 
the dwarf irregular galaxy Holmberg II. For the latter we assume 
an essentially flat metallicity gradient with a slope of only 1%, 
since dwarf galaxies typically have no metallicity gradients, but 
the slope cannot be exactly zer o if a dust-to-gas / metal s gradient 
is to be obtained from the model. Moustaka s et alj a2010f) also pro- 
vide gradients based on abundances derived using the P-mefhod 
(Pilvugin & Thuan 2005), which suggest significantly lower metal- 
licities and somewhat flatter gradients. Using the latter, the amount 
of metals are clearly inconsiste nt with the dust abundances derived 
by|M unoz-Mateos et al.l d2009l) . i.e., dust-to-metals ratios are above 
unity for a majority of data points. Therefore, we present only mod- 
els (see secti on|372) bas ed on m etallicity gradients derived using the 
iKobulnickv & Kewlevl ( 120040 calibration, which are also closer to 
the dust-to-gas gradients in terms of steepness (cf. the models with 
pure stellar dust production in Figure Q] which in several cases are 
relatively close to the dust-to-gas gradients). 

This inconsistency (that the dust-to-metals ratio exceeds 
unity), is indeed a problem, and the reader should bear in mind 
that we have not used these metallicity data because they are the 
most reliable, but rather because they are not obviously incon- 
sistent with the dust-to-gas ratios we compare with. As will be 
shown below, there are still a few data points at small galacto- 
centric distances which would correspond to dust-to-metals ratios 
ab ove unity also when using the m etallicity gradients derived using 
the lKobulnickv & Kewlevl j2004h calibration. Either the dust abun- 
dances should be about a factor of 2-3 lower, or the metal abun- 
dance should be increased by a similar factor. The latter is less 
likely, however, since we already use oxygen abundances which 
are at the high end of the possible range (set by the various mefh- 
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ods and calibrations for abundance determination) to estimate the 
metallicity along the disc and assume a relatively low oxygen frac- 
tion. Hence, it is more likely the dust abundance is overestimated 
in general, which is an interpretation that is supported by recen t 
results on the dust content of dwarf galaxies (Madde n et alj|201lh 
and which we will discuss in more detail in section [331 But for 
now, in order to avoid constructing strange inconsistent models, we 
have just scaled down the dust abundances for some galaxies such 
that the dust-to-metals ratio never exceeds 0.9 anywhere along the 
discs. 



3.2 Models and fitting 

We have fitted simple analytical models (see section \23\ 
to the dust-to-gas profile s derived from observations by 
M unoz- Mateos et all (2009). In doing so we have used the 
Levenberg-Markwardt scheme for ^-minimisation. More pre- 
cisely , we have used the IDL-routine package MPFIT ( Mar kwardt I 
2009) in combination with a numerical implementation (for IDL) of 
the Kummer-Tricomi functions (see Paper I). This setup has proven 
very stable for the cases we consider here (equations|10|andim). 



3.2.1 Input data 

For the model fits we have fixed the metal yield y z to the value 
obtained from a simple closed-box scenario 



Tz 



Z(R = 0.4R 25 ) 



(12) 



ln(l/ju) 

where \i is the overall gas mass fraction of the galaxy (see Table[U(. 
The metallicity at a galactocentric distance R = 0.4R 2 s is known 
to be a good proxy for the typical metallicity of a galaxy disc 
dGarnett et al.ll2002l) . which is why we define the effective metal 
yield as above. We have used the same gas mass fractions as in 
IPilvugin. Vflchez & Continil J2004I) . Note that the stellar dust yield 
ya, on the other hand, is treated as a free parameter. 

Furthermore, we use oxygen (O/H) gradients deriv ed using the 
strong-line calibration by Kobulnickv & Kewlev (2004) as input to 
our model (see lMoustakas et ai.ll2O10L for further details). We use 
the log-linear-fits and do not consider possible bends and features 
of the metallicity profiles as metallicities derived from H n regions 
can be quite uncertain. To obtain the total metal fraction (metallic- 
ity) we first convert the number a bundances of oxyge n into oxygen 
mass fractions using the relation ( Garnett et al.ll2002l) 



X = 12-, 
H 



(13) 



in which we then have implicitly assumed M g!is = 1.33Mh. Then, 
we assume that oxygen typically makes up about a third of all met- 
als (this is certain ly at the low end of the possible range, see e.g., 
iGarnett et af]|2002l. where 45-60% is the suggested value), hence, 



(14) 



The fraction of oxygen in the metals is of course not a universal 
constant and may vary from one environment to another. In partic- 
ular, the oxygen fraction is likely higher in low-metallicty environ- 
ments than it is in high-metallicity environments where low- and 
intermediate stars as well as SNe type la have contributed signifi- 
cantly to the metal enrichment of the ISM. 



3.2.2 Modelling results 

Allowing dust destruction (6 > 0) does not seem to improve the 
agreement between model and the dust-to-gas profile inferred from 
observations. For example, we considered the case where 6 is fixed 
to a specific value (we assumed 6 = 10, see section 1231 for a moti- 
vation of this value), but with e being an essentially free parameter 
constrained only by e > 6/Z max (where Z max is the maximum/central 
metallicity of the galaxy) in order to avoid the singularity that oc- 
curs at Z = 6/e (see Paper I). In such case the dust-to-gas/metals 
profiles cannot be reproduced for any of the galaxies in our sample. 
With 6 as a third free parameter (but still e > (5/Z max ), the best fit 
is typically obtained for 6 < 1. The only exceptions are NGC 2403 
and 3031 for which moderate dust destruction (5 ~ 4) leads to a 
slight improvement, albeit with a too high degree of freedom for 
the model fit to be better from a statistical point of view. This more 
or less rules out any significant net destruction of dust. Hence, we 
present here (in detail) only models without dust destruction due to 
SN shocks for the selected galaxies, since dust destruction hardly 
improves the models but rather just adds another free parameter 
(constraining more the two free parameters is also clearly not justi- 
fied from a statistical point of view, given the small number of data 
points we have at our disposal for several galaxies). 

Dust growth in the ISM, on the other hand acts as to improve 
the overall agreement between models and dust-to-gas/metals pro- 
files derived from observations. In fact, allowing dust growth (e > 
0, 6 = 0) improves the fits rather dramatically compared to the 
case of pure stellar dust production (e = S = 0). We find only one 
case (NGC 3031) where a really good fit cannot be obtained. This 
galaxy show a peculiar central depression in the dust-to-gas profile 
(see figureQ] middle panel on the second row from the top), which 
may not have anything to do with the physics of dust formation. 
A possible scenario is that a merger-induced starburst may have 
cleansed the central regions from dust, since NGC 303 1 appear in 
a relatively dense environment (it is for example interacting with 
NGC 3034 and 3077). For all other galaxies (including the outer 
part of NGC 3031) we obtain more or less good fits (see figure Q] 
and figure [2} . 

In most cases there certainly seem to be need for dust growth, 
or another kind of 'secondary' dust production, that depends on 
dust abundance and/or metallicity and thus steepens the dust-to- 
gas gradient and creates a dust-to-metals gradient. The very steep 
dust-to-metals gradients in the outer parts of NGC 2403 and 303 1, 
as well as that of NGC 3351, are particularly good examples. In 
most cases the required stellar dust yield y<j is well below 50% of 
the total metal yield y z according to the model fits with 6 = and 
6 being free a parameter (see Table [TJ. This is a very reasonable 
result, suggesting our model is sound. Models with 6 = 0, but with 
a significant e and relatively modest (or even insignificant) stellar 
dust production pr ovide very good fits in genera l to the dust-to-gas 
profiles derived bv lMunoz-Mateos et a 1. (2009). The cases requir- 
ing a high stellar dust yield have high dust-to-metals ratios over the 
whole disc (after down-scaling of the dust abundance), suggesting 
that not only may the dust abundances derived from observations 
for these galaxies be too high, but also the metallicities could be 
too low. 

In figure [3] we show e as function of the 5-band 
luminosity (data taken from IPilvugin & Ferrinil l2000l: 
IPilvugin. Vflchez & Continil l2004t) and the mean gas mass 
density, here estimated as 



<S g ) 



M„ 



2nR\; 



(15) 
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(data taken from |Pilvugin 1 Vflchez & Continil 120041; IWalter et al.l 
l2008l;lMunoz-Mateos et al.l2009l;lMoustakas et alj20ld) . formod- 
els where the fitting does not imply e -» (NGC 5055, 5194). 
The e-values correlate nicely with the B-band luminosity suggest- 
ing that dust-growth in the ISM is much more efficient in smaller 
and less luminous late-type galaxies. It is unclear why this correla- 
tion exists, but we believe there could be a connection to star forma- 
tion and the formation of cold molecular clouds. Low-luminosity 
disc galaxies tend to have less star formation which may create a 
suitable environment for dust growth in molecular clouds. If star 
formation is too intense, dust growth may simply be inhibited by 
the radiation field. Indirect support for this hypothesis is that pri- 
marily young stars (due to their emission at UV/blue wavelen gths) 
may be responsible for heating cold dust dGordon et al1l2004l) . The 
e-values does not at all seem to correlate with (X g ), however, which 
suggest the parameterisation we introduced in Paper I is very much 
reasonable. 



3.3 Comments on models of individual galaxies 

Here we summarise and comment on the modelling results of indi- 
vidual galaxies. 



3.3.1 NGC 628 

The dust-to-metals gradient is clearly negative. Assuming pure stel- 
lar dust production implies that the stellar dust yield ya must be 
equal to the total metal yield y z in order to match the dust abun- 
dance in general, which is obviously unrealistic (especially since 
the dust abundance has been scaled down by 62%). Adding signif- 
icant dust growth in the ISM provides a good model with y d being 
only about 10% of y z . A fairly strong indication of non-stellar dust 
production. 



3.3.2 NGC 925 

This galaxy also show a negative dust-to-metals gradient, but the 
slope is shallower. The dust-to-gas/metals gradient can be nicely 
modelled with dust growth. May be (or have been) interacting with 
a low-mass galax y, which could explain dynamical and morpholog- 
ical asymmetries iPis ano. Wilcots & Elm egreen 1 9981) . 



3.3.3 NGC 2403 

This is a somewhat peculiar case. The inner part of the dust-to- 
gas gradient has a slope which is very much consistent with the 
slope of the metallicity gradient, while the outer part is significantly 
steeper. Beyond R/R^s ~ 0.3 the observationally inferred dust-to- 
metals profile show a steep drop which is clearly inconsistent with 
the model with only stellar dust production. In that case, the model 
prediction falls outside the error bars for R/R25 > 0.7. 

Adding dust growth improves the agreement between model 
and observation tremendously. This is one of the strongest cases 
in support of grain growth being the dominant mechanism for dust 
production in the ISM. However, the very inefficient stellar dust 
production (yd only ~ 1% of yz) suggested by the model fit to the 
data is indicative of a possibly too low dust-to-metals ratio in the 
outer disc. The ratio is less than 1% in the outskirts of the disc and 
~ 50% in the inner disc (see figure[2]l. 



3.3.4 NGC 2841 

The dust-to-metals gradient is quite steep, but flattens out towards 
the centre of the galaxy. Both features can be perfectly reproduced 
by a model including dust growth and the required/implied stellar 
dust yield is just 7% of the metal yield, which makes this galaxy an- 
other strong case for dust growth - in particular since the error bars 
are quite small in this case. A model fit with only stellar dust pro- 
duction is clearly inconsistent with the dust-to-gas/metals gradient. 
The dust abundance was only scaled down by 5% in this case. 

3.3.5 NGC 3031 

NGC 303 1(M 81) is another peculiar case, where there is a clear 
dust depletion in the central parts. This is certainly not predicted by 
the model. The mid part of the disc has a dust-to-metals gradient 
which is nearly flat, but beyond R/R25 ~ 0.6 it steepens in a way 
similar to NGC 2403. It is worth noting that this galaxy is known to 
be in a relative ly dense enviro nment and interact with at least two 
other galaxies (Davidge 2008), which suggest that the central parts 
likely have formed by one or more mergers. A merger is a violent 
event, quite capable of destroying large amounts of dust, which is 
a possible explanation of the central depression in the dust-to-gas 
profile. SN-shock waves may not be efficient dust destroyers, al- 
though the environment in merger-driven nuclear starburst may be 
hostile enough to actually destroy dust grains and not just crush 
them into smaller grains. A model with significant dust growth in 
the ISM, but almost no stellar dust production, reproduces the slope 
of the outer dust-to-gas/metals profile nicely. Again, we argue this 
is strong evidence for grain growth in the ISM being the most im- 
portant dust formation mechanism. The dust abundance was scaled 
down by 25%. 

3.3.6 NGC 3198 

Once again the dust-to-metals gradient is clearly negative. Assum- 
ing pure stellar dust production is quite reasonable in this case (y<j 
being about 36% of y z ), although including dust growth enhances 
the agreement between model and observation significantly. Stellar 
dust production is still relatively important in this case, where y& is 
about 19% of y z . 

3.3.7 NGC 3351 

This case is interesting, since it shows the largest discrepancy be- 
tween the distributions of metals and dust. The metallicity gradi- 
ent is almost flat allover the disc. A model without dust growth 
is ruled out for two reasons: pure stellar dust production implies 
ya = yz, and the corresponding dust-to-gas/metals profile does not 
fit the data at all. Adding dust growth, however, we obtain a clearly 
better fit to the observations (see figure Q] and figure [2j- Yet an- 
other case in favour of somewhat lower dust abundances. To en- 
sure the dust-to-metals ratio never exceeds 0.9, the dust abundance 
had to be scaled down by a massive 85%, w hich indicate there 
might be a serious problem with the method (see lDraine & Lil2007l : 
iMunoz-Mateos et al.ll2009l) for deriving dust abundances from ob- 
servations. 

3.3.8 NGC 3521 

The model with pure stellar dust production fits quite nicely to the 
observational dust-do-gas/metals profile, albeit with a rather large 



© 201 1 RAS, MNRAS 000,[TJ-?? 



6 Mattsson & Andersen 

Table 1. Parameter values for the models. All models have 6 = 0. 



Galaxy 


Type 




e = 




£> 






Remark 






yz 


yi 


yilyz 


yi 


yilyz 


e 




NGC 628 


SAc 


0.0314 


0.0324 


1.00 


0.00322 


0.103 


120 




NGC 925 


SABd 


0.0156 


0.0029 


0.18 


0.00151 


0.097 


68.0 


Interactions with low-mass galaxy? 


NGC 2403 


SABcd 


0.0188 


0.0051 


0.27 


0.00000 


- 


457 




NGC 2841 


SA(r)b 


0.0292 


0.0107 


0.37 


0.00201 


0.069 


94.4 




ngc 303 i 


SAab 


0.0133 


0.0047 


0.35 


0001 8 


0.014 


216 


Interacting with NGC 1014 1077 


NGC 3198 


SBc 


0.0148 


0.0054 


0.36 


0.00287 


0.194 


39.6 




NGC 3351 


SBb 


0.0310 


0.0315 


1.00 


0.00000 




240 


Strong bar 


NGC 3521 


SAB be 


0.0200 


0.0150 


0.75 


0.01149 


0.575 


24.9 




NGC 3621 


SAd? 


0.0166 


0.0105 


0.63 


0.00376 


0.227 


177 


Pure disc/bulgeless galaxy 


NGC 4736 


SA(r)ab 


0.0119 


0.0069 


0.56 


0.00001 


0.001 


152 




NGC 5055 


SAbc 


0.0306 


0.0189 


0.62 


0.01825 


0.596 


0.00 


Interacting with UGC 8365? 


NGC 5194 


SAB be 


0.0247 


0.0217 


0.88 


0.02470 


1.000 


0.00 


Interacting with NGC 5195 


NGC 7331 


SAb 


0.0190 


0.0068 


0.36 


0.00586 


0.308 


9.94 


Retrograde bulge 


NGC 7793 


SAd 


0.0144 


0.0092 


0.64 


0.00284 


0.197 


140 




Holmberg II 


dlrr 


0.0086 


0.0005 


0.05 


0.00000 




14018 


No measured metallicity slope 



stellar dust yield (yd being about 75% of y z ) and with an exception 
for the inner most data point. With dust growth, the inner most point 
is still not reproduced, so the agreement with the observations is not 
clearly better (see figure Q] and figure [2}. The dust abundance was 
scaled down by 71%, which again provide evidence that there is a 
problem with how dust abundances are derived. 



3.3.9 NGC 3621 

The dust-to-metals gradient have clearly a negative slope, but the 
disagreement with a pure stellar dust production is moderate. Dust 
growth improves the overall agreement between model and obser- 
vation, and an essentially perfect fit is obtained for very reasonable 
parameter values. The model clearly reproduces the dust-to-gas as 
well as the dust-to-metals profile within the error bars. The dust 
abundance was scaled down by 35%. 

NGC 3621 is considered to be a pure disc (bulgeless) galaxy, 
which imply a quiescent formation history without major merg- 
ers and that central re gions are fairly u nevolved compared to other 
late-type galaxies dGliozzi et alj |2009). Hence, this galaxy may be 
a good test case for an idealised model such as the one we have 
applied here. 



3.3.10 NGC 4736 

This galaxy has a clear dust-to-metals gradient, which is nicely re- 
produced by a model with significant dust growth. Stellar dust pro- 
duction is essentially rejected by the fitting algorithm (y d ~ 10~ 5 ), 
which is unphysical, although it implies dust growth is needed and 
clearly dominating. A model fit with only stellar dust production 
is inconsistent with the dust-to-gas gradient, but almost within the 
error bars, which makes it marginally consistent with the dust-to- 
metals gradient since the error bars are larger for dust-to-metals 
ratios. However, a model with dust growth is clearly preferred. The 
dust abundance was scaled down by 54%, which yet again seems 
problematic. 



3.3.11 NGC 5055 

Here, the dust-to-metals gradient is almost exactly flat (see figure 
O. As a consequence the fitting routine reject the dust growth con- 
tribution (e — » 0). The required stellar dust yield vj is 62% of the 
metal yield y z , which is not a totally unreasonable figure. The inner 
most data point suggest a mild central depletion in the dust-to-gas 
profile. NGC 5055 is possibly interacting with UGC 8365. The dust 
abundance was scaled down by 14%. 

3.3.12 NGC 5194 

NGC 5194 is the only galaxy in the present study which deviates 
from the general trend that dust-to-metals gradients have negative 
slopes. In this galaxy the dust-to-gas profile is in principle flat over 
the whole disc, while the dust-to-metals gradient shows a mildly 
positive slope. Hence, NGC 5194 could seem better modelled with 
some moderate dust destruction, but dust destruction would require 
a higher stellar dust yield, and since the model with only stellar dust 
production suggest that 88% of the metals ejected from stars is in 
the form of dust, there is not much room for dust destruction. Dust 
growth means no improvement of the fit in this case. This galaxy 
show unrealistic dust-to-metals ratios at several locations along the 
disc, without down-scaling of the dust abundance. The dust abun- 
dance was scaled down by 60%, which was not enough to bring 
down the stellar dust yield to a reasonable number, however. This 
galaxy is somewhat unusual in that it appears to have an essentially 
flat metallicity gradient, something which is not expected for this 
morphological type (SABbc). The need for a 60% down-scaling 
of the dust abundance would not be if the metallicity of the cen- 
tral parts were as much higher as needed to create a more typical 
metallicity gradient. 

3.3.13 NGC 7331 

The dust-to-metals gradient is again nearly flat. Adding dust growth 
means only a minor improvement and the stellar dust yield required 
in the pure stellar dust model is very much reasonable (yd being 
about 36% of y z ). NGC 733 1 is therefore an example of that there is 
not always an obvious need for dust growth. A remarkable property 
of this galaxy is its retrograde bulge (counter-rotating relative to the 
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Figure 1. Dust-to-gas ratio as function of galactrocentric distance (circles) together with the best-fit models with stellar dust production only (dashed blue 
lines) and simple models including dust growth (solid red lines). The grey symbols in the background shows the original dust-to-gas ratio before lowering the 
dust abundance. 
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Figure 2. Dust-to-metal ratio as function of galactrocentric distance (circles) together with the best-fit models with stellar dust production only (dashed 
blue lines) and simple models including dust growth (solid red lines). The dotted (black) line shows the case of 50% metals locked-up in dust, which roughly 
corresponds to the dust-to-metals ratio in the Solar neighbourhood. The grey symbols in the background shows the original dust-to-metals ratio before lowering 
the dust abundance. The light-red/pink shaded regions correspond to dust-to-metals ratios above unity. 
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Figure 3. The correlation between the dust-growth parameter e (for models where the fitting does not imply e — > 0) and the B-band luminosity (left) and 
average gas mass density (right) of each galaxy. 



disc), a fact which is not likely to be connected to the distributions 
of metals and dust. 



3.3.14 NGC 7793 

This galaxy shows a moderately negative dust-to-metals gradient. 
A model with only stellar dust production works reasonably well, 
although it is not totally within the error bars. Adding dust growth 
yields an excellent model fit. The dust abundance was scaled down 
by 15%. 



3.3.15 Holmberg II 

This galaxy is the only one in our sample which is not really a disc 
galaxy. It is a dwarf irregular galaxy, which seems to have an essen- 
tially flat metallicity gradient and generally a low dust abundance, 
as is usually the case with galaxies of this morphological type. 
There is a dust-to-metals gradient, however, which can be nicely 
reproduced by a model with dust growth, provided the metallic- 
ity gradient is not exactly flat. Here we have (arbitrarily) assumed 
a 1% slope, which gives a reasonable result (although the high e- 
value makes Holmberg II deviate from the overall location of the 
data points in the right panel of figure[3]l. The required/implied stel- 
lar dust yield is approaching zero (y d — > 0), suggesting dust growth 
must be the dominating process for dust production. However, it 
should be noted that a model with only stellar dust production is 
in principle consistent (within the errors) of the dust-to-gas/metals 
gradients derived from observations. The stellar dust yield would 
in such case need to be only 5% of the metal yield. 



3.4 If dust is not destroyed - then what is the effect of shocks 
from SNe? 

Clearly, our results presented above suggest dust is not destroyed 
in any significant quantities due to SN shocks, as have been sug- 
gested by several studies and authors (see, e.g. graine & Salpeter 
19791; lMcKeelll989l; pHelens et alj[l994l: iNozawa & Kozasd 120061 
Jones, Tielens & Hollenbach] 1 19961: Ijonesl |2004 IJones & Nuthl 
201ll) .We propose dust reprocessing rather than destruction, which 



we motivate as follows. 

The first thing to consider is what we mean by 'dust destruc- 
tion'. Here, we prefer to look at dust destruction as a process that 
will atomise a dust grain/molecule. A SN-shock wave contains a 
vast amount of energy, which can certainly be harmful for dust 
grains, but it is quite possible that dust grains are shattered to small 
splinters, rather than completely evapourated. Small grains tend to 
survive better when exposed to e nergetic shock waves (see, e.g. 
IJones. T ielens &Hollenbachll996h . which implies there might be a 
lower size limit at which the shattering stops being effective. Such a 
scenario has the advantage of being favourble to dust growth in the 
ISM. Assuming spherical grains the total surface area A tol of all dust 
grains is related to the total volume V tol of dust as A tot 
where (a) is the the average grain radius. Thus, if the average 
grain size is decreases by a certain factor, the total surface area 
A tot increases by the same factor. The dust-growth timescale then 
decreases with a similar factor (see equation [6). Recent, more de- 
tail ed, analysis of th i s fact is presented by I Jones & Nuthl d201ll) 
and Jgirashita & Kuol fcOllh . In the latter study. iHirashita & Kuol 
d201 lb show that the time-scale of the grain growth is very sensitive 
to the grain-size distribution, because the grain growth is mainly 
regulated by the surface-to-volume ratio of grains, as mentioned 
above. 

This 'grain-shattering scenario' also helps to avoid a too large 
population of extremely large grains in the ISM, as recent de- 
velopments in modelling of dust growth and wind formation in 
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AGB-star atmospheres suggest dust grains may be r elatively large 
(micron-sized) already w hen they leave these stars dHofner 2008; 
iMattsson & Hofner 201 1). Attempts to solve the inverse-problem 
of finding the grain-size distribution from interstellar extinction 
measurements sugge st dust grains ar e typically not larger than 
~ lyum (see, e.g lClavton et al.l l2003). which would pretty much 
rule out the existence of extremely large grains in the ISM. AGB 
stars contribute a significant fraction of stellar dust, so the appar- 
ent absence of grains larger than ~ 1/jm must either suggest that 
grains are both grown and shattered to smaller fragments in the 
ISM, or the micron-sized grains predicted by dynamic AGB-star 
atmosphere modelling is a false prediction. 

Dust grains may also have their chemical-bond structure al- 
tered due t o external forc e s rathe r than being shattered or de s troyed 
(see, e.g. iTielens et al] 1 19871 ; Ijones, Dulev & Williamsl Il99d : 
IVollmer et alJl2007T) . For example, production of graphite-like car- 
bon dust and even nano-diamonds, is a likely outcome of shock- 
processing of interstellar dust. This type of dust re-processing is de- 
sirable since it is unlikely stars can produce the amount of graphite 
(or graphite-like) material needed to explain the 2175A-bump in, 



e.g., t h e M il ky W a y extinction cu rve (Mathis, Rumpl & Nordsieck 



19771; IPeil 1 19921 ; iKim, Martin & Hendrvl 1 1994 IClavton etal 



2003). Carbonaceous dus t may also produce hyd rocarbons when 



exposed to shock waves dTavlor & Williamsll 19931) . which would 
be desirable for sim ilar reasons, as hyd rocarbons may contribute to 
the 2175 A-feature ((Bl anco et al. 1991) as well as explain infrared 
features dKwok & Zhanghoilh . 

From a more general point of view, there is evidence for dust 
reprocessing in SN re mnants from diff e rences in dust species ra- 
tios. In a recent study, lAndersen et a found the ratio of 
very small grains to big grains to be higher than that found in the 
plane of the Milky Way in a study of molecular interacting SN rem- 
nants. The discrepancy is typically a factor of 23. They suggest 
that dust shattering is responsible for the relative over-abundance of 
small grains, which is consistent with the grain-shattering scenario 
we sugg est and in agreemen t with prediction from dust destruction 
models djones & Nuthl 1201 lh . In some cases small carbon grains 
are severely und erabundant, which is l ikely evidence for sputter- 
ing, according to lAndersen et alj d201ll) . 

3.5 Are dust abundances generally overestimated? 

As it was pointed out above, we believe the dust abundance may 
be ove restimated in a sev eral (if not all) cases . iMunoz-Mateos et al.l 
d2009h uses the model by|5 raine &Lild2007h which assumes a spe- 
cific composition of different dust species based on the inferred 
dust properties of the solar neighbourhood (except the variations of 
the PAH abundance) and a fixed grain-size distribution. We believe 
this to be a too simple model of the dust component in a galaxy, 
mainly because the dust properties in the local Galaxy most cer- 
tainly do not reflect the dust properties everywhere else and that 
the grain-size distribut ion most likely changes with environment 
fairashita&Kuoll201lh . 

We would also like to point out th at the dust-to - hydro gen ra- 
tio suggested for the local Galaxy by iDraine et alj f2007) , i.e., 
Ma/Ma = 0.010, is similar to the metals-to-hydrogen ratio for the 
solar neighbourhood (essentially the same as in the Sun), which 
implies a dust -to-metals rati o of ~ 0.6. Observationally, this ratio 
is ~ 0.3 - 0.5 j Whittedfl 99 lh which is also rough ly the ratio foun d 
in most local galaxies of similar type (see, e.g. Ilssa et al1ll990h . 
Hence, there is a possibility that the dust abundance could be over- 
estimated due to the model itself. Evidence in support of this idea 



has recently been presented by Mad den et al.1 J20T J), suggesting 
that the dust-to-gas ratios of dwarf galaxies are not in accord with 

the low metallicity of these systems. 

Furthermore, the estimates by IDraine &Lil d2007l) uses 

a total carbon abundance of 245 + 30 ppm relative to H 

I 1 J — — 3 1 1 

from Asplund et al. (2005). However, Landi, Feldman & Doschek 

(2007) and ICenteno & Socas-Navarrd d2008h obtain values of 
(0/H)q which are ~ 1.9 times la rger than the solar ox ygen value 
of (O/H)© = 457 + 56ppm from I Asplund et alj d2005l) . The solar 
carbon abundanc e might therefore b e significantly larger than the 
245ppm used bv lDraine&TH l feOOTl) . If the total carbon abundance 
is rather 350 pp m, the mass of carbon in dust w ould increase by 
75 % according to Shen. Draine & Johnson! ( f2008h . Note, however, 
that ICompiegne et al. 1( 1201 ]]) also obtain M^/Mh = 0.010 using 
a different model than that of IDraine & Lil d2007l) . which suggest 
a total carbon abundance of 269 + 33 ppm relative to H (from 
lAsplund et ai] |2009) where 26 % is gas phase carb on (to be com- 
pared to the 4% in the model by IDraine & Lill2007l) . 

Even if the dust abundances may need to be scaled down, we 
believe the errors are systematic in a way which will not affect the 
actual slope of the dust-to-gas profile very much. Hence, we con- 
sider the results we have obtained here being most likely qualita- 
tively correct, i.e., the dust-to-gas gradients are indeed steeper than 
the metallicity gradients, suggesting significant dust growth in the 
ISM. 



4 CONCLUSIONS 

Dust destruction by shocks from exploding SNe is believed to be a 
rather efficient process, evaporating dust grains such that large vol- 
umes of surrounding interstellar space is essentially cleared from 
dust. We challenge this viewpoint, and suggest that SN shocks in- 
duce reprocessing of interstellar dust rather than destroying it. We 
base our conclusion on the fact that the effect of dust destruction by 
SN-shock waves on the dust-to-metals gradients in late-type galax- 
ies is not compatible with these gradients being negative in general, 
which is what observations tell us at this point. SN-shocks may 
indeed still crush dust grains into small splinters, thus changing 
the grain-size distribution, but perhaps not disintegrate dust grains 
completely. Moreover, dust grains may change properties in other 
ways as well, e.g., through annealing or restructuring of the chemi- 
cal bonds. The latter may be a way to produce graphite-like dust in 
the ISM - a type of carbon dust which is unlikely to be formed in 
any larger quantities in stellar atmospheres. 

Two other results also stand out: (1) dust growth in the ISM 
seem to be the most important dust formation mechanism, and (2) 
conventional methods for estimating dust abundances may include 
systematic errors which are not fully under control. Assuming stel- 
lar dust production only, we found that dust-to-gas gradients are 
generally too steep to be reproduced by such a model. Dust growth 
is an efficient way to steepen the dust-to-gas gradient relative to the 
metallicity gradient, which makes it an inevitable ingredient in any 
attempt to model the distribution of dust along galactic discs. In 
many cases we found the dust growth to be the totally dominant 
dust production mechanism. However, the need for dust growth 
may be slightly exaggerated in our models, since there are cases 
suggesting the dust abundances derived from observations are too 
high, even after the dust abundances have been scaled down. In 
these cases the dust-to-metals ratio (as obtained from observations) 
exceeds unity without rescaling, which is clearly unphysical. If the 
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amount of dust is systematically overestimated, this may of course 
affect the need for dust growth. 

Regardless of whether the dust abundances should be lower 
or not, the metal content of the inner part of the galaxies seem 
dominated by dust in most cases. Even if the dust-to-metals ra- 
tios suggest the dust abundances may be overestimated in the cen- 
tral regions in particular, there is a general trend towards higher 
dust-to-metals at low galactocentric radii (a natural consequence of 
the dust-to-gas gradient being steeper than the metallicity gradient). 
This suggest the extinction properties may vary significantly over 
the disc of a late-type galaxy and the amount of extinction may 
not simply scale with metallicity. If so, corrections for extinction 
should not be made using a generic correction method applied to 
the whole disc. 
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